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Abstract 



Asymptotic expressions for Clebsch-Gordan coefficients are derived from 

an exact integral representation. Both the classically allowed and forbidden 

regions are analyzed. Higher-order approximations are calculated. These 

CS ' give, for example, six digit accuracy when the quantum numbers are in the 

hundreds. 
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I. INTRODUCTION 

This paper contains a detailed study of the asymptotics of Clebsch-Gordan coefficients 
and includes the derivation of new results. We use the term "Clebsch-Gordan coefficient" 
in its colloquial sense, i. e. the vector addition coefficients of SU{2). Thus our results also 
give the asymptotics of the 3j-symbols. We consider the case in which all of the quantum 
numbers get large together. What this means is multiplying all of the quantum numbers by 
a number and studying the asymptotic behavior of the Clebsch-Gordan coefficient as this 
multiplier gets large. Such a multiplier is often called 1/^, so that the limit of large quantum 
numbers is the limit of small h. 

The history of this subject dates back to the early days of quantum mechanics and the 
study of the classical limit of quantum mechanical quantities. Numerous papers have been 
written in this area. We summarize the literature briefly here. In 1959, Wigner [|I| discussed 
the physical interpretation and classical limits of Clebsch-Gordan coefficients. He described 
a certain average behavior, and did not analyze the oscillatory nature of the Clebsch-Gordan 
coefficients. There are references in this work to Edmonds and Brussaard and Tolhoek 
0]. In 1968, Ponzano and Regge [Q presented asymptotic expressions that included the 
oscillations. Their work included an interpretation of certain angles that occur in their 
results and in ours. Additionally, they discussed the allowed and forbidden regions. However 
their derivation is, in their words, "rather heuristic." It was borne out in their comparisons 
with the exact values. William Miller p[ derived similar expressions using semiclassical 
methods in 1974, but did not treat the forbidden region. Another work that relates to the 
present paper is that of Srinivasa Rao and V. Rajeswari 0. It contains exact expressions for 
Clebsch-Gordan coefficients and their relationship to certain hypergeometric series. There 
is more information in the work of Biedenharn and Louck . 

In this paper, we start by deriving an exact integral representation for the Clebsch- 
Gordan coefficients. Then the methods of stationary phase are used to approximate this 
integral. The allowed and forbidden regions are treated separately, and the resulting ex- 
pressions are related to the literature. These methods are then used to derive higher-order 
results, that is, the next order in an expansion in h. These formulas are accurate to five or 
six digits when the quantum numbers are in the hundreds. 

Possible applications of this work include high-angular momentum calculations and the- 
oretical investigations which contain sums over large numbers of Clebsch-Gordan coefficients 



II. EXACT EXPRESSIONS FOR THE CLEBSCH-GORDAN COEFFICIENT 

Our starting point is an exact expression for the Clebsch-Gordan (vector-addition) coef- 
ficient, due to Wigner (see, for example, Eq. (3.6.11) of Ref. 0), 

{jimiJ2m2\jm) = [{2j + l)(ji+J2-j)Kji-J2+j)K-Ji+J2+j)V(ji+J2+j + l)!]^^^ 

X {(ji+miy.(ji-miy.{J2+m2)\{J2-m2)\{j +m)\{j -m)\Y^'^ 

.. y^ ^~^y . (2.1) 

' ^ rl I Ti -I- T,^ — i — r ill Ti — m, — ^\Ui^-i-mn — r ill t — in -I- m, -I- 7- ill t — i-, — m^ -I- r il 



A factor of Sm,mi+m2 has been omitted; throughout this paper we will assume that m is equal 
to mi + 7712- Also, unless otherwise specified, sums over an index are sums over all integers. 
It will turn out, though, that the summand is nonzero for only finitely many values of the 
index. 

We begin by deriving the following exact expression for the Clebsch-Gordan coefficient. 



(ji mi J2 m2 I j m) 



-lY+^N. 



]i rax J2 ni2 J m 



(ji - mi)\ {J2 - m2)\ 



duj \dt 



j2-in2 



[{t - iy+J2-j'i (t - uy^+^'~^ {u - iy+^^-^^] \^__ 



■0,t=0 ' 



(2.2) 



where Nj^rnij2m2jm is defined to be 



TV, 



Jimij2m2jm 



(2j + l){ji+mi)\{ji-mi)\{J2+m2)\{J2-m2y.{j + m)\{j-m) 
(ii +J2 +3 + 1)1 iji +J2 - j) ! iji -J2+j) ! (-J1 + J2 +j) ! 



1/2 



(2.3) 



Because the quantity being differentiated in Eq. ( |2.2|) is a polynomial in the variables u and 
t, the operation of differentiating this quantity and then evaluating the result at u = and 
t = simply selects a particular coefficient in the polynomial. Thus, Eq. ( p.2| ) expresses the 
Clebsch-Gordan coefficient as a certain coefficient in a polynomial that can be written in 



closed form. This equation can be derived from results in the literature |jTO|], but we give here 
an independent derivation of Eq. (|2.2| ) from Eq. ( p.l|) to verify that all of the conventions 
involved are consistent. 

In order to prove Eq. (|2.2|) , we start by finding the coefficient of u^^~"^'^ P^'"^^ in the 
polynomial (t — l)J+^2-ji (j^ _ y^yi+ji-j (^^ _ \y+ji-j2_ This is equal to the coefficient of 



u 



n-m 



1 in the polynomial that is given by {u — 1) 



j+ji-32 



times the tt-dependent coefficient of 



p2 m,2 jj^ ^]-^g polynomial (t — 1)-J+-J2 h ^f _ y^yi+J2 i_ Using the Binomial Theorem, we get 



(t - 1)^+J'2-Ji = Y^ 



j + J2 - ji 
k 



t^(-l) 



j+J2-ji-k 



(2.4) 



and 



(*-")-*-- =E(-"^rO 



The coefficient of P'^ "^^ in the product of these is 



J +32-31 
k 



-u 



\jl+J2-j- 



iy+j2-ji-k I 31+32-3 
32 -m2-k 



-u 



\3i+32-j-{J2-m2~'k) 



(2.5) 



= (iy2+m2 ^h-3+m2 Sp ( 3+32- Jl \ / Jl + J2 " 3 

^^ \ k J \j2-rn2- k 

As explained above, we need to multiply this polynomial by 



u 



(2.6) 



j + jl - J2 



u _ iy+h-j2 = J2( I U' i-iy^''-^'-' . (2.7) 



and find the coefficient of u^'^ ™^ . Tlie result is 

(_l)i2+m2 Sp f J+J2- Jl \ f Jl+ 32 -3 \ (3+31- 32 \ .^y,-j,+m+k 

^ \^ k J \J2 — 1^2 — kjyj — m — kj 

= f^iy+m \^(-iY ( J + -^2 - ii ^ ( Jl +32-3\ ( 3 + 3i - 32 \ /2 8) 

^ ^ ^^ ' \32 + m2-z J \ z J \ji-mi-z J ' ^ ' ^ 

where we have redefined the index of summation according to k = j — ji — m2 + -z in the 

final line of this equation and made use of the identity I t I = I ^ ) • Since j ' + m is 

always an integer, (— l)2(i+m) jg equal to one, and we have shown that the right-hand side 
of Eq. (p^ ) is equal to 

N,,m.,,^,,my{-iy{'V'~'' 1 (^'^^^-^\ (3+31-32 \ 

nmxnmijm /_^\ ) y^^j^^^_ ^ j \ z ) \3i-rni- z ) 

{-ly Nj^ mi jamaim (j + h " Jl)! (jl + 32 " j)! (j + Jl " J2)! 



E 



(j - Jl -m2 + zy. (J2 + m2- zy (jl + J2 -j - zy.zl (jl -nil- zy (j - J2 + mi + zy 

(2.9) 



This is the same as the right-hand side of Eq. (p.l|) and completes the proof of Eq. (^.21). 
An alternative proof begins by introducing a factor of x^ into the sum in Eq. ( p.l|) and 
deriving a third-order differential equation for the resulting function of x. This differential 
equation can be solved using hypergeometric functions, and the result eventually leads to 
the expression shown in Eq. ( p.2| ). 

Equation (|2.2|) can be used to obtain an exact expression for the Clebsch-Gordan coef- 
ficient as an integral. One uses the orthogonality of the functions exp(m^) on the interval 
[— TT, vr] to select the desired coefficients in the polynomials. Thus, we substitute exp(i6') for 
t and exp(i0) for u in the polynomial (t - 1)^+^2-^1 (t - uy^+^^'^ (u - ly+i^-^^ in Eq. (^, 
multiply by exp[— ^(ji — mi)(j) — i{32 ~ "^2)^]; and integrate the two variable from — vr to vr. 
The resulting expression for the Clebsch-Gordan coefficient is 

(jl mi ]2 ms I j m) = (-l)-'+™ A^^^ 






TV rn 



— 7r ^ —IT 



(2.10) 

This may be rewritten using the definition of the sin function, whereupon it becomes natural 
to redefine the angles by a factor of two. The resulting form is 



(ji mi J2 m2 \jm) = (-1)^+'" {2iy+^'+^^ tt'^ A^, 



Jimij2m2jm 



/n/2 i-k/2 
/ e^'"''^^^'"'^^ sin-'+^'^-^'^ e sm^'+^^-^{e - 0) sin^+^'i-^'^ rf^ rf0 . (2.11) 

-7r/2 i-7r/2 



It is this integral expression for the Clebsch-Gordan coefficient that we use in the following 
sections to derive formulas for the asymptotic behavior of these coefficients. 

It is also possible to express the Clebsch-Gordan coefficient as a coefficient of a term 
in a polynomial in one variable, and thus as a one-dimensional integral. Equation ( p.2|) 
shows how the Clebsch-Gordan coefficient is related to the coefficient of u-^i"™-! p2-'^2 j^ ^]^q 
polynomial {t — ly+^^-ji (f _ ^^ii+j2-j (^^ _ iy+ji-j2_ xhis is the same as the coefficient of 
uji-m^ ^^My2-m2 jn ^i^g polynomial {u^^ - iy+^^~J^ {u^^ - uy'+^^~^ {u - 1)^+^1-^2 (that is, 
t has been replaced by u'^^) for sufficiently large integers M. This can be seen as follows. 
We start by imagining the polynomial (t — l)J+J2-ii ^^ _ y^yi+j2-j ^y^ _ iy+ji-j2 expanded 
out into a sum of monomials. If t is replaced by u^^ , each of the monomials is now just 
a coefficient times a power of u. We do not want any of these terms to have the same 
power of u, otherwise they would combine and the coefficients would change. Thus, we 
look at the original polynomial (t — l)J+J2-ji ^f _ y_yi+j2-j (^^ _ -iy+ji-j2 g^j^j g^gj^ what the 
highest power of u is. This is (ji + J2 — j) + (j + ji — J2) = 2ji. We therefore select 
M to be 2ji + 1. The result is that the coefficient of ^i-Ji^™! tJ2~™-2 in the polynomial 
{t - l)J+J2-ii (t _ u)Ji+J2-i ^y _ iy+ji-j2 ig the same as the coefficient of u^i-"^i+(2ii+i)(i2-m2) 
in the polynomial (u^-^i+^ — l)-^+-^2-ii (^^2ji+i_y^ji+i2-i (^y^_iy+ji-J2_ ^^ j^g^y ^^j.^^ g^j^ overall 
factor of u^^~^^^~^ , so this coefficient is the same as the coefficient of -uJ-^+^ii 02-^2) j^^ ^^le 
polynomial (m^-^^^^ — l)J+J2-ii ^y^2ji _ ]^yi+j2-j ^y^ _ \y+h-32 _ "YYv^ resulting expression for 
the Clebsch-Gordan coefficient as a coefficient in a polynomial in one variable is 



(ji "^1 h m2 \jm) = (-1)^'+"^ A^, 



Ji mij2m2j m 



[j - m + 2ji{J2 - 7712) 



(1 \ ]-m+Z]\\.]2-'m.2) 
— r(u2^i+^ - i)^^+J2-ii (y2n _ iy,+32-j (y _ ly+h~n^ 1 _ . (2.12) 

duj '"-" 

As above, the selection of the coefficient in the polynomial can also be carried out with an 
integral. 



{jimiJ2m2\jm) = {-ly^"" Nj^rniJ2m2Jm 



1 

2^ 



X / g-*b-™+2ii(i2-m2)]0 /gj(2ji+i)?!> _ -^y+j2-ji (^i'2ji(t> _ ]_~\ii+i2-j (^i4> _ -iy+h-J2 ^A /'2 i^^ 

This may be rewritten as 

C2^)i+ii+i2 
(ji mi j2 "12 I j m) = ( - 1 )■' +™ Nj^ rm n ^2 j m — 



^i{2nm2+m)^ si^J+J2-n [(^-^ + 1/2)0] sin^'i+J2-J(ji 0) sinJ'+J'i-J2(0/2) d0 , (2.14) 

■n 

and this may be simplified to the form 

(jimij2m2 I jm) = (-l)^'+'"%^,,,^,,^7r-i 2^+^'^+^'^ (2.15) 

X / cos[(2jim2 + m)(f) + ^{j + ji + J2)] 

X sin^'+^'2~^'^[(ji + 1/2)0] sin^'^^^-^iji 0) sin^'+^'^'-''2(0/2) # . 

Although this is a one-dimensional integral (as opposed to the two-dimensional integral 
presented above), it seems to be not as useful for the study of asymptotics because of the 
presence of the magnetic quantum numbers in the argument of the cosine function. 



III. STATIONARY-PHASE APPROXIMATION OF INTEGRAL EXPRESSION 
FOR THE CLEBSCH-GORDAN COEFFICIENT 

In order to carry out a stationary-phase approximation of the integral expression for 
the Clebsch-Gordan coefficient presented in the previous section, we begin by writing the 
expression in Eq. ( p.ll|) in the form 



{jim,j2m2\jm) = {-iy-^"'{2zy+^^+^'n-'N,,^,,,^,,^ T T e'^'''^U9d(P . (3.1) 

J-Tr/2 J -tt/2 

where the function g{6, 0) is defined to be 

g{e, 0) = 2zmi0 + 2im2e + (j + J2 - Ji) ln(sin 6) (3.2) 

+ (ji + 32 - 3) ln[sin(6' - 0)] + (j + ji - J2) ln(sin 0) . 

Note that g{6, 0) has singularities where it goes to —00, but the integral is still well-defined 
because the integrand is exp(f7). 

To find the stationary-phase points (as explained in Appendix A), we must first compute 
the first derivatives of the function g. 

da 

-^ = 2im2 + (j + 32 - 31) cot e + (ji + 32 - 3) coi{9 - 0) , (3.3) 

da 

-^ = 2imi - (ji + j2 - 3) cot(^ - 0) + (j + ji - J2) cot . (3.4) 

Setting these first derivatives equal to zero results in a system of two equations in two 
variables. The identity 

^fa A\ l + cot^cot0 
cot^ - = — — —- 3.5 

cot — cot fe* 

may be used to transform this system to an equivalent system. 



2im + {] + J2 - Ji) cot e+{i+ ii - J2) cot = (3.6) 

/. . -N /I / ■ ■ ., 1 + cot ^ cote/) „ ,„ „, 

2^m2 + (j + J2 - Ji) cot^ + {j,+j^-j)^— -f = 0. (3.7 

cot (p — cot tl 

In order to be clear on phase conventions, choices of signs and branch cuts we write out 
the steps involved in solving this system of two equations for cot 6 and cot 0. We start by 
multiplying the second equation by (j + ji — J2)(cot — cot 9) and substituting in first one: 

[2im2 + U + ]2 - Ji) cot 9]{-[2im + (j + j2 - Ji) cot 9] - (j + ji - J2) cot 9} 

+ (Ji + J2 - ]){{] + Ji - 32) - cot 9[2im + (j + j2 - Ji) cot ^]} = 0. (3.8) 

This is a quadratic equation in cot 9. 
cot^ ^[(j + j2 _ j,)(_2j) - {31+32 -3){3 +32 - Ji)] + cot 9[{] + j2 - Ji)(-2im) + 2zm2(-2j) 

+ (ii + J2 - j)(-2«"^)] + 2zm2(-2im) + (ji + j2 - j)(j + Ji - J2) = . (3.9) 

Simplifying this results in 

-COt2^(j+j2-Jl)(jl+j2+j) 

- 4i cot 9{i2m + ms j) + 4m2m + (ji + J2 - j) (j + j'l - J2) = . (3. 10) 

The two solutions for the quantities cot 9 and cot are (the upper choice of sign is one 
solution and the lower choice of sign is the other). 

^n -2i{J2m + m2J) T P 

cot 9 



ijl + 32 + j){j + J2 - jl) 



, , -2i{jim + mij) ± p 

COt(p=— : —-. : — , I'J-ll) 

[31 + 32 + 3){3 + 31 -32) 

where /3 is defined to be 

j3 = ^ 4:17111712 f - Ammiil - 4:mm2Jf + (ji + J2 - j)(j + h - 3i)U + Ji - J2)(ji + J2 + j) ■ 

(3.12) 

In this equation, we use the usual choice of branch cut for the square-root function: if 
the argument is negative, then the result is a positive number times the imaginary unit. As 
discussed in Sec. [Ill A|, the quantity P is real for classically allowed sets of quantum numbers. 



and it is pure imaginary for classically forbidden sets of quantum numbers. It should be 
noted that this is the same definition for the symbol P as in Ref. |^. 

The stationary-phase approximation of the integral /T^ /g Jli/2 e^^^''^^ d9 dcj) that appears 
in the expression for the Clebsch-Gordan coefficient in Eq. ( p.l| ) is given by a sum of terms 
of the form 



^"^ :e^(^'^), (3.13) 



det 



suinmed over stationary-phase points. The branch cut for the square root function is just 
below the negative imaginary axis, as is usual. The symbol Q7gjx2 denotes the 2x2 Hessian 
matrix of second-order derivatives of the function g[6, (p), whose entries are given by 

U + h - ji) csc^ 9 - (ji + J2 - j) csc^(6' - (j)) , 



89^ 

090(1) 

Ocj)'^ 



(jl+j2-j)cSc2(^-0), 

-(ji + J2 - j) csc'(^ - 0) - (j + Jl - J2) CSC^ . (3.14) 



Using the identity csc^ 9 = 1 + cot^ 9 these quantities can be expressed in terms of the 
cotangents in Eq. ( 3.11|) without addressing the issue of branch cuts of the arc-cotangent 



function. The value of csc^(^ — 0) can be determined from the quantities in Eq. ( p.ll|) using 



the identity sin(6' — 0) = sin 6* cos — sin cos 6* = sin 6* sin0(cot — cot 6*). The determinant 
becomes 

O'^g 
d^t 0//1 A^9 = U + 32 - ji) csc^ 9{j + ji - ja) csc^ 



0(9,, 



+ (ji + J2 - j) csc^{9 - 0) [(j + J2 - ji) csc^ 9 + {j + ji - ja) csc^ 0] , 

= (1 + COt^ 9){1 + COt^ 0){(j + J2 - J,){j + Ji - J2) 

+ J;^^'^:;'j/gp [(j +J2 - Ji)(l + cot^^) + (j +Ji - J2)(l + cot^0)]}. (3.15) 

In this form the determinant is expressed entirely in term of the cotangents of 9 and 0. The 
quantity e^^^'^^ can also be expressed in this way. Choices of branch cuts are not necessary 
when expressing sin 9 and sin in terms of the cotangents because only even powers of the 
sine- functions appear. Using the identity 

e'^ = cos9 + ism9 = sin 9{cot 9 + i) , (3. 16) 

we obtain for the factor e^^^'^' in Eq. (|3.13|) 

= {i + cot 0)2™i(i + cot ^)2™2 sin^-+i2-ii+2m2 Q sinii+i2-i(^ _ 0) ^:^^3+3i-J2+2m, ^ 
= {i + cot 0)2™i {i + cot ^)2™2 sin2^2+2m2 Q ^^^^ ^ _ p^^ Qyi+j2-j sjj^2ii+2mi ^ 

= {i + cot0)2™Hi + cot^)'™^ (1 + cot2^)-^'2-'"^(cot0 - cot^)^'^+^'2-^' (1 + cotV)-^'i-'"i 
- (^ + -^^r-^-^ (. + cot.)-- (eot0-cot.)^-^--- (3.17) 



{-i + cot 0)Ji+'^i {-i + cot 6')J2+'"2 

Using this equation and Eq. ( p.l5| ), all of the quantities in the expression in Eq. ( |3.13| ) can be 
expressed in term of the cotangents of 9 and 0, given in Eq. (|3.11| ). It should be noted that 



all of the exponents in Eq. (|3.17|) are integers, so choices of branch cuts are not necessary. 



A. Allowed region 

It is useful to introduce the concepts of a triangle-allowed region and a classically allowed 
region of the space of values for the quantum numbers. We define the triangle- allowed region 
to be the set of quantum numbers for which ji, J2 and j satisfy the triangle inequalities and 
for which the inequalities \m\ < j and {|wij| < ji, i = 1,2} hold. The Clebsch-Gordan 
coefficient is zero outside of this region, so it is only within this region that asymptotic 
expressions are desired. The triangle-allowed region is divided into a classically allowed 
region and a classically forbidden region. As is usual, we call these the allowed and forbidden 
regions for brevity. The allowed region is defined to be the set of quantum numbers for 
which it is possible to define j-vectors in a three-dimensional space in such a way that their 
lengths are equal to the j-values and their 2;-components are equal to the m-values (and, 
of course, such that j = ji + J2)- An example of such a construction for a set of allowed 
quantum numbers is shown in Fig. 1. It follows from the definition that the allowed region 
is contained in the triangle-allowed region. Examples of classically forbidden points are 
easily found in extreme cases, such as mi = ji. In this case, there is only one classically 
allowed value for m2 (assuming a set of triangle-allowed j-values have been given), because 
the ji-vector must point in the 2;-direction, and thus the j-triangle lies in a vertical plane. 

The allowed region is the same as the region in which the three A-values defined in 
Eq. ( |3.25| ) satisfy the triangle inequalities. This is because the A-values are the lengths 



of the projections of the j-vectors into the xy-plane. If the A-values satisfy the triangle 
inequalities, then it is possible to draw a triangle in the xy-p\a.ne with sides equal to the 
A-values. From this, one can construct the j-vectors by simply including the rn-values as 
z-components. Conversely, if the j-vectors can be constructed, then their projections into 
the xy-p\a.ne form a triangle (with the tail of J2 at the tip of ji), and the A-values satisfy the 
triangle inequalities. 

It is explained later in this paper that the A-values satisfy the triangle inequalities if and 
only if the quantity (— Ai + A2 + A3)(Ai — A2 + A3)(Ai + A2 — A3) is nonnegative (that is, it is 
not possible for two of the factors to be negative). This observation together with the fact 
that the quantity jS defined in Eq. (|3.12|) may be written as 



/3 = v/(Ai + A2 + A3)(-Ai + A2 + A3)(Ai - A2 + A3)(Ai + A2 - A3) (3.18) 

leads us to the result that the sign of /3^ distinguishes the allowed and forbidden regions: it 
is positive in the allowed region, and it is negative in the forbidden region. In the allowed 
region, (3 is four times the area of the triangle whose sides are the A-values. This triangle 
is the projection of the j-triangle into the xy-plane (see Fig. 1). From Eq. (|3.12|) it is 
apparent that for fixed values of the j quantum numbers, /?^ is a quadratic polynomial in 
the m quantum numbers. Thus, in the {1711,1712) plane the boundary between the allowed 
and forbidden regions is an ellipse. This is shown shown in Fig. 2 for one choice of values 
for ji, J2 and j. The boundary of the triangle-allowed region is the irregular hexagon. The 
forbidden region is composed of six subregions. The points that separate them are indicated 
in Fig. 2. These are the points where the ellipse that separates the allowed and forbidden 
regions is tangent to the hexagon that defines the triangle-allowed region. The coordinates 
of these points can be calculated from the expression for (3 and the equations for the straight- 
line sections of the boundary of the triangle-allowed region. The resulting coordinates of 

9 



these points are indicated in the figure. 

The calculations involved in the stationary-phase approximation of the integral expres- 
sion for the Clebsch-Gordan coefficient are different in the allowed and forbidden regions. We 
will treat the allowed region first. The sum over stationary-phase points for the case where 
the set of quantum numbers is in the allowed region is a sum over both of the two solutions 
for the cotangents of 9 and (f) given in Eq. (|3.11|) . This is analogous to the behavior demon- 



strated in Appendix A, and it is also the same as in the calculation of the stationary-phase 
approximation of the Airy integral, which is the canonical example of a stationary-phase 
calculation. In the case of the Airy integral, there are allowed and forbidden regions in 
position space, and in the allowed region the contour of integration is deformed to run over 
both of the stationary-phase points. We define ce and c^ to be the first solution for the 
cotangents in Eq. ( p.ll| ). 



Ce 



-2i{J2m + m2J) - (3 

ij +J2-jl){jl+J2+j) 



-2i{jim + mij) + (3 
[31 + 32 + 3){j + 3i ~ 32) 

Because of the form of the two solutions given in Eq. (|3.11|) , the second solution is obtained 
from this one by multiplying by minus one and complex conjugating. Note this is only 
valid in the allowed region, where the quantity (3 is real. The first term in the sum over 
stationary-phase points is given by plugging the expression for det g.g f^ , given in Eq. ( p. 151 ), 



and the expression for e^^^''^\ given in Eq. ( p.l7| ), into the quantity in Eq. ( |3.13|) , using cq 



and C0 for the cotangents. The second term in the sum over stationary-phase points is the 
same, except — c^ and — c^ are used for the cotangents. The result for the determinant in 
the second term is obtained by simply complex conjugating the first value, since all of the 
cotangents in this expression are squared. As for the e^^^''^^ factor in the second term, we 
start by considering the expression for this factor in the first term: 



The complex conjugate of this is 



icl - cl) 



*\jl+J2-j 



This shows that the e^^^'^^ factor in the second term [which appears after the (— l)Ji+J2+i 
in the last line of Eq. ( p.21|) ] is (— l)Ji+-J2+i times the complex conjugate of the e^^^'^' factor 



in the first term. Thus the second term is (— l)-?i+-J2+j times the complex conjugate of 
the first term. It is therefore convenient to obtain the sum over stationary-phase points 
by taking i-?i+.?2+i times the first term, adding the complex conjugate of this product, and 
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then dividing by i^i+-?2+i^ Using the fact that the real part of a quantity x is given by 
^[x] = {x + x*)/2, our stationary-phase approximation for the integral expression for the 
Clebsch-Gordan coefficient in Eq. (p.l|) can be written as 



(ji mi j2 m2 \jm)^ (_l)i+- (2ty+^^+^' n~^ N,, „, 



Jimij2m2jm 



iJi+J2+j 



^ 



271 in+h+j 



M9,<t>) 



det 






/_i\j+mr,i+ii+J2+2 -1 A^. . . so 



j^jl+J2+j 



o9{6,4') 



det 






(3.22) 



where the quantities det 



wg fp and e^^^'''^'' are obtained from Eq. (|3.15| ) and Eq. ( p.lTj ) using 
the ce and c^ given in Eq. (|3.19|) . 



Although the expression in Eq. ( |3.22| ) gives a value that is a real number, it involves 
intermediate quantities that are complex. It is possible to transform this expression so that 
only real quantities are involved. This transformation is very lengthy, and it is not practical 
to describe it in detail here. Instead, we present an expression that is exactly equal to the 
expression in Eq. (|3.22| ) in the allowed region. This equality can be verified most convincingly 
by substituting numerical values into the expressions and evaluating the results to high 
numerical precision (much higher than the level at which discrepancies would occur if order 
h terms were dropped). A brief description of the transformation is the following. Every 
complex quantity x + iy that occurs in Eq. (|3.22|) is written as the product of a modulus and 
a phase, a/x^ + y"^ exp[ztan~-'^(y/x)], where care must be taken that correct branches are 
used for each x + iy, that is, one must examine the quantities x and y to determine the range 
of phase factors {x + iy)/ yx^^-y^ that can occur in the allowed region, and make branch 
choices accordingly. At some stages in the calculation, large polynomials are involved, and 
computer-aided symbol manipulation becomes useful in working with these. Our result may 
be put in the form 



( jl mi J2 7712 \J7n) ^ 21 j^ rni J2 ma j m 

where x is defined to be 




cos 



vr 



X + ^ - vr(j + 1) 



(3.23) 



X 



Jl + - ) cos" 



+ 1 J2 + 2 I COS 



+ 1 J + 2 ) ^°s" 



{-7n){jl + jl - j^) - 7n2{jl + f - jl) 

a\i 
miif + Jl - jf) - (-m)(j| + jf - f) 

aX2 
m2ijf + f - jI) - mi{f + jl - ]l 



aA.q 



Till cos 



-|- 7712 COS 



-1 



-^1 + -^3 ~ -^2 
2A1A3 

-^3 + -^2 ~ -^1 
2A2A3 



(3.24) 
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and 



a= VU + h + 32) (-J + 31 + 32) (j - Ji + J2) (j + Ji - J2) 
Ai = A A? - m? z = l,2 

A3 



a/ j^ — r/i^ 



(3.25) 



(The quantity a is four times the area of the j-triangle shown in Fig. 1.) Note that the cos ^ 
functions in Eq. (|3.24|) are the usual principal branch, whose range is the interval from zero 



to TT. The quantity /jimijamajm 
-'ii mi J2 m2Jm = 



is defined to be 



(J + 1/2)(J+Ji+J2) 



3(3+31+32 + 1) 
fiji + mi) fiji - mi) f{J2 + ma) /(J2 - ^2) /(j + m) f{j - m) 
fiji + 32 + 3) fiJi + 32 - 3) fiJi - 32 + 3) f{-3i + 32 + 3) 

where the function / is defined to be 



(3.26) 



f{n) 



n\ 



V2 



imn"- e 



(3.27) 



that is, f{n) is the square root of the ratio of n\ to the Stirling approximation of n\. Note 
that for large n, f{n) approaches one. Thus, for large quantum numbers, Ij^mij2m2jra 
approaches one. It differs from unity by a correction that is order h, as can be deduced 
from the discussion of the Stirling approximation in Appendix A. As mentioned above, we 
present our approximation in the form given in Eq. ( 3.23|) so that the exact equality of this 
expression and the complex expression given in Eq. ( p.22| ) can be verified numerically. The 
may be dropped without reducing the quality of the approximation, that 



factor /, 



J\ m\j2'm2j m 



is, the ratio of our approximation to the exact value differs from unity by a quantity that is 
order h. Thus, we may write our approximation in the form 



( ji rni J2 m2 \ j m) 




X 



vrlJ + - 



(3.28) 



Ponzano and Regge [^ give a geometrical interpretation of the five angles that occur in the 
expression for x in Eq. ( |3.24D . An equation similar to Eq. ( |3.28| ) also appears in Ref. 0, but 
the (j + 1/2) factors in x are included at the end of the calculation to improve the accuracy, 
and the 7r(j + 1) in Eq. ( p.23| ) is missing so that the formula gives the wrong sign for even 
j'-values and does not give the right magnitude for half-integer j values. 



B. Forbidden region 

In the forbidden region, only one of the stationary-phase points is used in the approxi- 
mation. This is analogous to the situation in the Airy-function problem mentioned above, 
where there are dominant and subdominant branches, and in the forbidden region only the 
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subdominant branch exists. Similarly, the model problem in Appendix A shows how for the 
case oi m > n, two stationary-phase points are used, while for the case m < n only one 
stationary-phase point is involved. The choice of which of the two roots in Eq. ( p.ll| ) is to be 
used for our approximation of the Clebsch-Gordan coefficient is indicated in Table I. Given 
the m-values of a point in the forbidden region in Fig. 2, it is inconvenient to determine 
which subregion it is in by using nested if-then statements, because the relative ordering of, 
say, the m2-coordinates of the points on the boundaries between the forbidden subregions 
changes as the j-values are changed. A much simpler way to determine which branch to 
use is to find the sign of a certain polynomial which we describe here. As can be seen from 
Table I, the choice of branch alternates as one goes around the diagram in Fig. 2. Thus we 
use the sign of the product of three expressions that fiip signs in the right way. Given the 
coordinates of one of the boundary points in the (mi,m2)-plane, a vector perpendicular to 
it can be constructed by exchanging the coordinates and changing the sign of one of them. 
The dot-product of this vector and (mi, m2) is a function on the (mi, m2)-plane that changes 
sign at the boundary between the two subregions in question. Thus we are led to consider 
the sign of the function 

[(mi,m2)-(-2j|,/-j2_j2)] 

X [(mi, m2) ■ i-f + Ji - ji, f + Ji - 32 
X [irm,m2)-{-f+jl+ji,2jl 



If this quantity is positive (negative), then the upper (lower) choice of root in Eq. ( |3.11D is 



used. Once the cotangents of the angles at the stationary-phase point are determined, the 
approximation of the Clebsch-Gordan coefficient can be evaluated from the expression 

{Jimu2m2\jm) ^ {-lY^"' {2zy^''+^'^ n~' N,,^,,,m,jm , ^^ • (3.29) 



All of the quantities needed to evaluate this expression were expressed in terms of the 
cotangents of the angles in Eqs. ( ^.15| ) and ( |3.17| ). It may be noted that in the forbidden 
region the cotangents become pure imaginary, as can be seen from Eq. ( p.ll|) . This behavior 
is similar to that in the model problem in Appendix A, where the angle suddenly jumps 
in terms of its real part (but the analogy is not perfect because in the model problem the 
cotangent is pure imaginary in both the region m > n and the region m < n). 

Evaluating Eq. (|3.29| ) results in a real value, although complex numbers are involved at 
intermediate steps. As in the case of our analysis in the allowed region, the expression may 
be transformed to a form that involves only operations with real numbers. This can be done 
in a way that parallels the previous calculation, with hyperbolic functions playing the role 
of trigonometric functions. The transformation involves choices of branch cuts and depends 
on which of the six subregions of the forbidden region one is working in. Thus there are 
six different all-real expressions for the forbidden region. In the interest of brevity, we will 
present only one of these here. In subregion VI, the expression in Eq. ( |3.29|) is exactly equal 
to 
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.]^y'2+m2 2 / 



31 mi J2m2j m 



vr 



exp(-X 



{vi)-\ 



(3.30) 



where y'-^*-' is defined to be 



X 



(vi) 



ji + - ) cosh ^ 



- ( J2 + 2 ) cosh ^ 



J + - ) cosh ^ 



"-"^(j? + Jl - f) - m2{jf + f - 


-j1) 


aAi 




"-"^i(j' + ji - fi) - m{ji + jf - 


-3') 


0X2 




-m2{jl + j2 - j|) + ^i(j' + Jl 


-jD 



aAq 



— -m cosh ^ 

— ?Ti2 cosh~^ 



-^1 + -^3 ~ K 
2A1A3 

2A2A1 



(3.31) 



This all-real expression was derived by a very lengthy calculation, as in the case of the 
analysis in the allowed region. Again, an exact equality such as the one above can be 
checked easily by substituting in test numbers and evaluating to sufficient precision. As 
before, to actually use the approximation, one would drop the factor of Ijimij2m2jm since 
it can be approximated by unity, to the order that we are working in this section. All-real 
expressions for the other subregions of the forbidden region can most easily be obtained by 
using the symmetries of the Clebsch-Gordan coefficients to related the expressions for the 
different subregions. If one prefers not to work with six different expressions for the forbidden 
region, one can use the polynomial discussed above to select the required stationary-phase 
point and then plug this into the approximation given in Eq. ( ^.29| ). This requires operations 
with complex numbers, but is easier to implement in a computer program. Alternatively, 
to obtain an approximate value for the Clebsch-Gordan coefficient for a given point in the 
forbidden region, one could work with only one all-real expression for a particular forbidden 
subregion and use the symmetries of the Clebsch-Gordan coefficients to map the given point 
to a point that is within the subregion for which the expression is valid. 

It is interesting to compare the all-real expressions obtained in the allowed region, 
Eq. ( p.23| ), and in the forbidden region, Eq. (|3.3CI| ). They are similar in form, but the 
behavior is oscillatory in the allowed region and exponentially decaying in the forbidden re- 
gion. This is the behavior expected in quantum mechanical problems that have an allowed 
region and a forbidden region. 

In the forbidden region, writing the approximation in an all-real form is illuminating 
because it makes it apparent that sign functions exist. We call the factor (— l)J2+™-2 jj^ 
Eq. (13.301 ) a sign function. The remaining factors in that equation are all positive, so the 
sign function gives the sign of the result. However, since the result is an approximation of 
a Clebsch-Gordan coefficient, the sign function also gives the sign of the Clebsch-Gordan 
coefficient, at least in the asymptotic regime. Thus, the sign functions are actually properties 
of the Clebsch-Gordan coefficients themselves, for a given choice of phase conventions. We 
are using the conventions defined by Eq. (p.lj) . The existence of sign functions was not clear 
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from Eq. ( p.l|) , which was our starting point. The sign functions for each of the six forbidden 
subregions are given in Table I. 

The angle Xi given in Eq. ( |3.24| ), that appears in our approximation in the allowed region 
can be rewritten in several different ways. The reason is that the angles that multiply the 
m's in the equation for x ^'^^ two of the interior angles in the triangle formed by the three 
A- values. If we call these angles ai, 02 and 0:3 (where ttj is the angle opposite the side of 
length Aj), then we have 

iTLi + m2 = rn , 
«! + ^2 + «3 = TT . (3.32) 

Thus, the vectors (mi, m2, — m) and (ai, 0^2, as — vr) are both perpendicular to (1, 1, 1), and 
their cross-product is parallel to (1, 1, 1). Each of the three components of their cross product 
are thus equal, and each one could be used as part of x in the allowed region, 

miQ;2 — m2ai = ^-2(0:3 — vr) + ma2 = —mai — mi{as — tt) . (3.33) 

In the forbidden region there is no such flexibility in how to write the corresponding terms, 
because there do not exist three angles corresponding to ai, 0:2 and 03. This is because it 
is not possible to form a triangle using the three A-values. Three different quantities like 

can be written down by cyclically permuting the indices, but only two of these can be used as 
arguments of the cosh" function in an all-real expression. This can be seen in the following 
way. The A's are nonnegative and if three nonnegative numbers fail to satisfy the triangle 
inequalities, exactly one triangle inequality is violated. [Proof: Let \max be the largest value. 
Amid be the middle value, and Xmin be the smallest. Then — A^m + ^mid + ^max > and 
Amm- Amid + Amax > 0, SO wc must havc Amm + Amid - Amax < 0.] Now we consider rewriting 
the expression 

Aj + A^ - Aj ^ ^ ^ (A3-Ai)^-A^ ^ ^ _ (Ai + A2-A3)(A2 + A2-Ai) ^^ ^^^ 

2A3A1 2A3A1 2A3A1 

This shows that of the three permutations of the expression in Eq. ( |3.34| ), exactly two will 
be greater than unity. It is these two that must be used as arguments of the cosh~^ function 
in an all-real expression. Thus there is no flexibility in ways to write the m-terms in x &s 
in the allowed region. Throughout each one of the six subregions of the forbidden region, a 
single triangle inequality for the A's is violated. It is not possible that one triangle inequality 
is violated in one part of a subregion and another triangle inequality is violated in another 
part of the same subregion because at the boundary between these two parts (3 would be 
zero, as can be seen from Eq. (|3.18| ). But (3'^ is a quadratic polynomial in uii and m2 [see 



Eq. ( |3.12| )], the zero-contour of which is the ellipse in Fig. 2, so it is not possible for it to 
be zero along a curve in the forbidden region. The A that is largest in each subregion is 
indicated in Table I. The forms of all-real expressions in each of the forbidden subregions 
will reflect the fact that in each one of the subregions one of the A's is larger than the sum 
of the other two. 
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It remains to discuss the case of points that are on the boundary between the allowed 
and forbidden regions. The quantity /3 in Eq. ( p. 12 ) is zero on this boundary, and since (3 



is invariant under the full 72-element symmetry group of the 3-j symbol [|T0| , the Clebsch- 
Gordan coefficient cannot be approximated on the boundary with the formulas presented 
in this paper. The reason is that (3 occurs in the denominator in Eqs. ( |3.28| ) and ( p. 301 ). 



Since /3^ is a homogeneous polynomial in the quantum numbers, it will be zero for sets of 
quantum numbers equal to any multiple of a set of quantum numbers for which f3 is zero. 
The behavior of the Clebsch-Gordan coefficients in the direction transverse to the boundary 
should be similar to that of the Airy function (see Ref. ||^). 

The invariance of (3 under the 72-element symmetry group of the 3-j symbols may be 
shown as follows. We begin by constructing the 3x3 Regge array of linear combinations 
of quantum numbers, given in Ref. [|Ty] . For any integer n, we define the polynomial pn to 



be the sum of the n-th powers of the nine elements of this matrix. These polynomials are 
invariant under the symmetry group, because if two sets of quantum numbers are related 
by a Regge symmetry, we can construct the 3x3 Regge array for each set and compute Pn. 
The results are the same because of the commutativity of addition. It is possible to write 
/3^ in terms of the p„. 

/?' = {pt - Qp2 pI - 27pl + 108p4)/324 . (3.36) 

The coefficients in this equation may be simplified slightly by using the relation pi = 3(ji + 
J2 + j). This equation proves the invariance of (3 under the symmetry group. 
An example of quantum numbers for which (3 is zero is 

{Jumij2,m2j,m) = (3,-2,6,4,7,2). (3.37) 

This point is not on the edge of the triangle-allowed region. Points for which (3 is zero and 
which are on the edge of the triangle-allowed region are easier to find. For example, one can 
choose (ji,mi) = {J2,fn2) = |(j, m). For such a point, ji + J2 — j is zero. 

The reason we are unable to approximate the Clebsch-Gordan coefficient for cases in 
which f3 is zero is that the determinant of the 2x2 matrix of second derivatives of g{6, 0) 
is zero at the stationary-phase points. This can be shown by plugging the solutions for 
the cotangents of 6 and at a stationary-phase point into Eq. ( 3.15| ) for the determinant; 



the result has an overall factor of (3 after being simplified [see Eq. ([B7|)]. This determinant 
appears in the denominator of Eq. ( |3.13|) , so our method cannot be applied. Note that when 



(3 is zero, the two solutions for the cotangents at the stationary-phase point are the same 
[see Eq. (|3.11| )]. Also, it should be noted that if any of the m- values has its absolute value 



close to the corresponding j, then the corresponding A will be small [see Eq. (|3.25| )], and the 
area of the A-triangle will be small. Thus, f3 will be small, and the set of quantum numbers 
is close to the boundary. In contrast to this, there is no difficulty with the approximation 
if the m values are close to zero. These considerations are mirrored in the approximation 
(using Stirling's formula) of Nj-^rnij2m2jm, defined in Eq. (|2.3|) ; no factorials of m- values 
appear, only factorials of j — m, j + m, etc. 
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IV. HIGHER-ORDER APPROXIMATION 

The methods used in the previous sections can be extended to higher order. In this 
section, we derive the next correction to the previous results. The approximation that is 
obtained in this way gives results that are accurate to six digits, for example, when the 
quantum numbers are in the hundreds. 

Let (6*0, 0o) be a stationary-phase point, i.e. a point at which gf = gf = 0. We write the 
Taylor expansion of the function g{9, (p) about the point [Oq, (J)q) as a sum of homogeneous 
polynomials, 

g{9o + x,(j)o + y) = go + 92 + 93 + 9i + ■■■, (4.1) 



where 



9o = 9{do,<Po) , 

92 = 9ee^'^l'^ + 9o4.^y + 9^yi'^ , 

93 = 9eoeX^/Q + 9eo^^^y/^ + 9e^^^y^/^ + £/.^^<pZ/V6 , 
^4 = 9eeeeXV'^'^ + 9eoe^^^y/^ + 9ee^^X^y^/'^ + 9e^^.Xy^/Q + 9^.^^y^ 1'^^ , (4.2) 



where, for example, 9^^^ is defined to be d^g/dO'^dcp at the stationary-phase point. 

To obtain the next higher stationary-phase approximation for the Clebsch-Gordan coef 



ficient, we terminate the series in Eq. ([4.1|) at the fourth-order term. The reason for this is 
explained below. Thus, the approximation of the function 9 has derivatives at the stationary- 
phase point {6q, 0o) that agree with those of 9 through fourth order. Our approximation of 
the integrand exp((yf) is 



exp[5f(6'o + X, 00 + y)] ~ exp(5fo) exp(5f2) exp(5f3) exp(5f4) 

= exp(f7o) exp(^2)(l + 93 + 91/^ + ...){! + 94 + 91/'^ + 



(4.3) 



When this is multiplied out, each of the terms may be integrated over the xy-plane in closed 
form. We are interested in the asymptotic behavior of the resulting terms. The question is 
how the terms behave when all of the quantum numbers {ji,Tni,J2,m2,j,m) are multiplied 
by the same factor (such a factor is called 1/h, as explained in the introduction). The 
stationary-phase point (^oi^o) is independent of the factor, i.e. (6*0,00) is order h^, as can 
be seen from Eq. ( p.ll|) . The second derivatives of 9 at (6*0, 0o) are order 1/h, as can be seen 



from Eq. (|3.14|) , so to see how the integral of exp(5f2) depends on h, we define new variables 
of integration to be the old variables times hr^'"^ . From this it follows that the integral of 
exp(5f2) is order h. By the same reasoning, the integral of a homogeneous quartic polynomial 
times exp(5f2) is order h^, and the integral of a homogeneous sixth-order polynomial times 
exp(5f2) is order h^. The polynomials 94, and 9I have coefficients that are order 1/h and 
l//i^, respectively, so the integrals of these times exp((72) are both order /i^. This is one 
order of h smaller than the integral of exp((72)- The integral of a homogeneous polynomial 
of odd degree times exp(5f2) vanishes due to antisymmetry. Thus, the next higher order 
approximation of the integral of exp(5f) is obtained by integrating 

exp(^o)exp(^2)(l + ^4 + ^3/2)- (4-4) 
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Terms coming from g^, etc, contribute at higher orders. 

To find the ratio of the integral of g4exp{g2) to the integral of exp{g2) the following 
integrals are necessary. 



r-OO /"OO 



h = 


/ 


/ exp{g2)dxdy, 




/ — OO .J 


'-OO 




j-CO 


/•OO 


i2 = 


/ 


/ x^ exp{g2) dx dy , 




/— OO .; 


l-oo 




j-CO 


/•OO 


h = 


/ 


/ x^yexp{g2)dxdy, 




/— OO .; 


'-00 




j-CO 


/•OO 


ii = 


/ 


/ x'^y'^ exp{g2) dx dy . 


1 


/~oo ^ 


'-00 

/ • 1 • / • ml 



(4.5) 

We will need the ratios i2/ii, is/ii and i4,/ii. The integrals are tabulated and these 
ratios can be worked out without the use of any information about relationships between 
the various derivatives of the function g at the stationary-phase point. The results are 

^2/^l 

"''"' ^ (9 9 -9'Y ■ ^ ^^ 

It may be noted that g^^ g — g^ is the determinant of the 2x2 matrix of second partial 
derivatives of the function g. 

The ratio, which we denote by ^4, of the integral of g4exp{g2) to the integral of exp(5f2) 
works out to be 



{9eo 


9^4 


. g2 )2 


- 


-^9e4 


9 ,j„f, 


i9go 


9^4 ~ 


.gl )2 

^04,' 


2^ + 9 


ee 9ij,(f, 



9eeee 9^^ ^9eg9tj> 9e4, y,/,^ ~r '^9eeg>4>\ 9g^ ' 9ee ",^0/ ^i/9000 9ee i/e<p ' 9444,4 9ffg 

8{9ee944-9gy 



(4.7) 



Next, we move on to the g^ term in Eq. ([4.4|) . To find the ratio of the integral of 
^glexp{g2) to the integral of exp{g2) the following integrals are necessary. 



H = 




/•OO 

/ x^ exp{g2) dx dy , 




/-OO J 


'-00 




/•OO 


/•OO 


ie = 


/ 


/ x^yexp{g2)dxdy, 




/-OO J 


'-00 




/•OO 


/•OO 


h = 


/ 


/ x^y^ exp{g2) dx dy 




/-OO J 


'-00 




/•OO 


/•OO 


is = 


/ 


/ x^y^exp{g2)dxdy 




/-OO J 


'-00 



(4.8) 

We will need the ratios i^/ii, ie/ii, ij/h and is/h- As in the case of the 5^4 calculation, the 
integrals are tabulated and these ratios can be worked out without the use of any information 



about relationships between the various derivatives of the function g at the stationary-phase 
point. The results are 



ySee 5^0<^ -'e</)'' 



"''"" ^ (9 9 -g'Y ' 

. ,. ^9o,{'^9l+39oe9j 

"''"' ^ {9 9 -9'f ■ ^ ^^ 

The ratio, which we denote by 5q, of the integral of \g1 exp((72) to the integral of exp(f72) 
works out to be 

•^6 = [^9e^{^9ge 9^^ + '^9lM9eoe 9^^^ + He^ 9e^^) 



- H9oe 9^^ + %L) (2^*89 9e^^ 9^4, + 2^^^^ 9^00 9 eg + ^9i^ 9 4^ + 3^L 9, 

,2 I „ „ 2 



ee^ci><t> ' Ho^/y HBBe iiB<t>tt, ^<t>ci> ' "ii<f>it,<f, ti,f,ee nee ' ■^ee<i,^4><t> -^4,(1,0^98} 

~r '-'^9gs\9e0e 9 004, 9j.j, ' 94,4,4, 94,40 9a 



^000 H004, ^4,4 ' H 444 H 440 H gg / 

'}444 9gg)\l[^^\9g0 944 iJg^j 



H9L 9l + 9L 9t)]/m900 944 - 9lf] ■ (4.10) 



The only remaining matter is to find the necessary values of the higher derivatives of the 
function g at the stationary-phase point. The definition of 5^, given in Eq. ( |3.3| ) can be used 
to find its second derivatives, given in Eq. (|3.14| ), its third derivatives, 

9^9 .,, , , , ^„„ 2 ,,„„,,, ,0.. , . .N„.„2 



d9^ 
d^g 

d^g 



2(j + 32 - Ji) esc' ^ cot ^ + 2(ji + j2 - 3) csc'(^ - 0) cot( 
-2(ji + 32 - 3) csc\9 - <P) cot(^ - 0) , 
2 2(ji + 32 - 3) csc\e - 0) cot(^ - 0) , 



dddct) 

^ = 2(ji + 32 - 3) csc2(0 - 0) cot(0 -e)+ 2(j + ji - 32) csc2 cot , (4.11) 

and its fourth derivatives, 



cfg 



-2(j + 32 - ji) esc' e (3 cot' e + 1 



- 2(ji + J2 - j) csc'(e - 0)[3cot'(^ - 0) + 1], 
^^ 2(ji + 32 - 3) csc'(0 - 0) [3 cot'(^ - 0) + 1] 



dO'^dcp 



0^9 



-2(ji + 32 - 3) csc'(^ - 0) [3 cot'(^ - 0) + 1] 



■ = 2(ji + 32 - 3) csc\e - 0) [3 cot'(0 - 0) + 1] 



i" ri/_- I _■ _-\ ,2/n J,\ [Q „^-(-2/ 
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= -2(ji + J2 - j) csc2(^ - 0) [3 cot2(e - 0) + 1] 

- 2(j + ji - j2) esc' (3 cot' + 1) . (4.12) 

Given the values of the cotangents of 6 and at a stationary-phase point, these derivatives 
can be evaluated without having to find the angles, i. e. without having to make any choices 
of branch cuts. One way to do this is to use Eq. ( p.5| ) to evaluate cot (6' — 0), and the identity 
esc' 9 = 1 + cot' 9 to evaluate the squared cosecants. 

The relationship between the Clebsch-Gordan coefficient and the integral of exp{g) is 
given in Eq. ( pj.ll ). Combining this with our higher-order approximation for the integral 
results in the following higher-order approximation for the Clebsch-Gordan coefficient. Each 
stationary-phase point contributes 






(_1)^-+™ {2t)^+n+n ^-2 N,,^,,,^,,^ , {I + 5^ + 5,) . (4.13) 



9g0 9^,4, 9g^ 



As explained in the previous sections, for points in the allowed region the sum over 
stationary-phase points is a sum over both of the solutions for the cotangents of 9 and 
(f) given in Eq. ( |3.11|) , and for points in the forbidden region only one of these solutions con- 



tributes. Given values for cot 9 and cot 0, Eqs. ( p.l4|) , ( [4.11| ) and (|4.12|) are used to evaluate 



the higher derivatives of the function g at the stationary-phase point. Then Eqs. ( [4.7|) and 
( [4.10| ) are used to obtain 64, and Sq. As shown in Eq. ( |3.17| ), the quantity e^*^^"'''^''^ can also 



be evaluated using the values of the cotangents of 9 and 0. The quantity Nj-^rnij2m2jm can 
be approximated to sufficient accuracy using the next correction to Stirling's approximation 
for the factorials. 

Finally, we present some numerical examples. We begin with the allowed region. 

For (ji,mi,J2,^2,j,^) = (200,100,300,150,400,250), the values are 

exact = 0.0703499 , 
approx = 0.0703496 . (4.14) 

For {ji,mi,J2,m2,j,m) = (200, 100, 300 + 1/2, 150 + 1/2,400 + 1/2, 250 + 1/2), the values 
are 

exact = 0.0730636 , 
approx = 0.0730633 . (4.15) 

In the forbidden region, the Clebsch-Gordan coefficients are much smaller. The following 
examples are from subregion I. 

For (ji,mi,J2,"^2,j,"^) = (200,150,300,-250,400,-100), the values are 

exact = 3.08961 x 10"^^ , 
approx = 3.08958 x 10"^^ . (4.16) 

For iJi,mi,J2,m2,j,m) = (200,150,300+1/2,-250 + 1/2,400 + 1/2,-100 + 1/2), the 
values are 
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exact = 5.32718 x 10"^^ 

-19 



approx = 5.32712 x lO"''^ . (4.17) 

Further examples of results from the higher-order approximation are discussed in Ap- 
pendix B. 



V. CONCLUSION 

The methods presented in this paper provide simple formulas for calculating first-order 
approximations to Clebsch-Gordan coefficients in the allowed region and in all of the for- 
bidden subregions. Additionally, a higher-order approximation is derived, although the 
expressions are more complicated. We do not know if the quantity ^4 + Sq in Eq. ( [4. 131) can 
be simplified when expressed in terms of the quantum numbers (see Appendix B for a special 
case). It appears to be complicated, as is often the case for higher-order approximations. 
The geometrical structure is not as clear. 

Our higher-order approximation provides the only known way to compute certain digits 
of some Clebsch-Gordan coefficients. By this we mean that given any computer, we can 
always find quantum numbers large enough so that the exact calculation is not feasible. 
The beginning digits may be calculated using first-order approximations; the higher-order 
approximation makes it possible to compute further digits. 

The methods of this paper could also be used to derive asymptotic expressions for the 
6j-symbols, etc. The starting point would again be an exact expression for the quantity 
of interest. One would then have to construct a polynomial with the property that the 
coefficient of one of its terms is this exact expression. Then an integral expression would 
be obtained, and finally this integral would be approximated using the stationary-phase 
method. 

As mentioned in the introduction, this work could have applications in high-angular mo- 
mentum calculations and theoretical investigations which contain sums over large numbers 
of Clebsch-Gordan coefficients [|,3 . 

Our analysis in the forbidden region led us to the realization that simple sign functions 
exist there that give the sign of the exact Clebsch-Gordan coefficients. These are summarized 
in Table I. 

A subject for future work is the approximation of Clebsch-Gordan coefficients and 6j- 
symbols near the boundary between the allowed and forbidden regions. Ponzano and Regge 
^ have conjectured and supplied numerical evidence for a typical Airy-function caustic 
behavior. Also, of course, it should be possible to extend the present calculations to even 
higher orders. 

APPENDIX A: A ONE-DIMENSIONAL EXAMPLE 

In this appendix we consider a one- dimensional example of an integral that gives a Fourier 
coefficient of a function which is an integer power of a fixed function. We are interested in 
the asymptotics of the result for large values of the two integers involved. 

We define the function F{m, n) for positive integers m and n by 
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I-k/2 

F{m, n) = cos" x e*™"^ dx . 

J-n/2 

It is possible to evaluate this integral exactly in closed form: 

F{m,n) = 



(Al) 



n 



n — 771 even 



{n — m)/2 

(_l)(n+l— )/2 2n+2 ^, [(m+n+l)/2]! (n^-n-l)! ^^ 

\ ' (m+n+lj! [(m— n— 1)/2J! ' 

2„+2^, [(n+l-H/2]![(n+l+H/2]! ^ _ ^ ^^^ > ^ 
(n+1— m)! (n+l+m)! ' 



(A2) 



In deriving these results, one uses the definition of the beta function, B{z + l,w + 1) = 
Jq t^(l — t)'^ dt, and the relation between the beta function and the gamma function, 
B{z,w) = T{z)r{w)/T{z + w). One also uses the results that for integers k >0, 



T{k + 



V^{2k)\ 

22k f.] 



r{-k + -) = {-!) 



.2^^^k\ 



(A3) 



1. Asymptotics of the exact expressions 

In order to compare the stationary-phase approximations derived in the following subsec- 
tion with the exact value of F{m,n), we will use Stirling's approximation for the factorials 
in the exact expressions in Eq. ( |A2[ ). The accuracy to which we will work is that the ratio of 
the exact value to the approximation should go to unity as n and m go to infinity, holding 
the ratio of n to m fixed. The difference between the logarithm of the exact expression and 
the logarithm of the approximation thus goes to zero as the two integers get large (the errors 
are of order 1/n). 

Stirling's approximation, through order unity (for the logarithms), is 



x\ 



V27rxx^e ^ 



(A4) 



The next correction to this is a multiplicative factor of e^'^^'^^\ Thus, the ratio of x\ to the 
approximation given in Eq. ( [A4| ) approaches unity as x goes to infinity. 
Our approximation of the exact expression for F{m, n) works out to be 



l_rn \ m/2 



1+- 



l-(f)^ 



-(n+l)/2 



n > m 



F{m,n) Ri < 



n — meven,n<m (A5) 



1 \(n+l-m)/2 2 [^ ( 



'-1^' 



' ™_i\ "*/2 



t+1 



rn\ " 
n ) 



-(n+l)/2 



n 



m odd, n < m 
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In deriving this result, we have used the fact that the inequahty n > m imphes n — m » 1. 
This is true because we are holding the ratio of the two integers fixed while letting them 
become large. In other words, errors of order 1/n are the same order as errors of order 
l/(n — m). Thus the Stirling approximation is used for quantities such as {n — m)\. Similar 
remarks apply to the inequality n < m. 

In Eq. (|A5|) , only positive quantities are raised to powers that could be non-integer. Thus 
there are no phase ambiguities. If one is sloppy about phases, the last expression appears 
to be the same as the first, differing only be a factor of two. The origin of this factor of 
two has a simple interpretation in the stationary-phase approximation, described in the next 
subsection. 

It is remarkable that the first expression (for the case n > m, n — m even) and the third 
expression (for the case n > m, n — m odd) in Eq. ( [A2| ) have the same asymptotics to the 
order at which we are working. A calculation is involved in showing this. The result that 
comes from applying the Stirling approximation to the third expression is 




m/2 



m 



n 



-(n+l)/2 



To the accuracy to which we are working, this turns out to be the same as the first expression 
in Eq. (|A5|) , although some work is required to show this. 



as 



2. Stationary-phase approximation 

To do a stationary-phase approximation for the function F{m,n), we write the function 

F{m,n) = / e"^(^)rfa;, (A6) 

where the function g is defined by 

g{z) =lncosz + i — z. (A7) 

n 

With the usual choice of branch cut for the logarithm function, the function g{z) is analytic 
everywhere in the complex plane except for vertical lines that intersect the real axis at odd 
multiples of tt, and at the intervals on the real axis where cos 2; is nonpositive. The identity 

cos(z -f iy) = cos x cosh y — i sin x sinh y (A8) 

is useful in showing this. Knowledge of the region of analyticity of g allows us to deform 
the contour of integration in Eq. ( |A6| ) without changing the value of the integral. We would 
like to deform the contour so that the phase of the integrand e""^^^' is constant. To do this, 
we need to know the imaginary part of g{z). With the help of Eq. { \KEi ) we find that this is 

'^[g{x + iy)] = — tan^"'^(tana; tanhj/) H x . (A9) 



n 
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If a contour is selected in such a way that this function is a constant, then g{x + iy) will equal 
an imaginary constant plus a real- valued function along the contour. The integrand e"^*^^-* 
will then equal a constant phase factor raised to the n-th power times a fixed real-valued 
function raised to the n-th power. This fixed real-valued function may be approximated by 
a Gaussian, and the integral may then be evaluated. Stationary-phase points Zg satisfy the 
condition 



9[Zs 



0, 



(AlO) 



This is equivalent to the condition 



tan2;„ 



.m 



n 



We note that the value of g" at a stationary-phase point is 



g {zs) = - sec Zs = ( — ) - 1 



(All) 



(A12) 



It is necessary to distinguish two cases, the case m < n and the case m> n (recall that 
m and n are both positive by assumption). We first consider the case m < n. In this case, 
it follows from a study of Eq. (|A9|) that a constant phase contour exists that connects the 
endpoints of the integral and passes through the stationary-phase point 



1 -1 "^ 

Zs = I tanh — 

n 

in a direction that is parallel to the real axis. The integrand is approximated by 

ng{zs)+ng"{zs){z-ZsY/2 



(A13) 



and the result for the integral is 



27r 



-ng"{zs) 



oiT-gizs) 




27r 
n 



1 m 

n_ 



m/2 



-o^ 



-1 -(n+l)/2 



(AM) 



which agrees with the result in Eq. 

We now move on to the case m > n. In this case, no single contour exists with the 
properties that it connect the endpoints of the integral and that the quantity in Eq. (|A9|) 
be constant. Instead, we choose a contour consisting of three straight-line pieces. The first 
part, Ci, is defined to start at — 7r/2 and go vertically upwards to —'7i/2 + iY, Y being a large 
positive real number. The second part, C2, is defined to go from —tt/2 + iV to tt/2 + iY, 
and the third part, C3, goes straight down to the 7r/2 endpoint of the integral. The parts 
Ci and C3 contain stationary-phase points, which we call Zs- and Zs+, and which are given 
by 



Zs± = ± — h t tanh — 
2 m 



(A15) 



For the integral along Ci the integrand is approximated by 
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^ng{zs^)+ng"{zs-)(z-Zs-)'^ /2 



We parametrize the curve Cihy z = Zg- + it, where t is a real parameter. Then dz is i dt 
and the resulting approximation for the integral along Ci is 



2n 



^ng(zs-) 



+ng"{zs-) 
Similarly, the approximation for the integral along C3 is 



(A16) 



27r 



— t 



+ng"{zs 



^ng{zs+) 



(A17) 



Because of the condition m > n the integral along the curve C2 goes to zero as Y goes 
to infinity. The approximation for F{m,n) is thus the sum of the expressions given in 
Eqs. ( [A16| ) and ([A17|) . If ?ti — n is odd the sum is zero, in agreement with the exact value. If 
m — n is even, the sum agrees with the approximation of the exact result, given in Eq. ( |A5D . 
Thus we see that depending on the ratio m/n, different numbers of stationary-phase 
points must be considered due to fundamental changes in the form of the stationary-phase 
contours as m/n goes from one side of the critical value of 1 to the other side. On either 
side of the critical value of 1, a stationary-phase approximation is possible. The behavior 
near the critical value is discussed briefly in the main part of this paper. 



APPENDIX B: THE CASE OF VANISHING MAGNETIC QUANTUM 

NUMBERS 



Clebsch-Gordan coefficients for the case of vanishing magnetic quantum numbers (mi = 
m2 = m = 0) are of interest in atomic and nuclear physics. Many of the expressions derived 
in this paper simplify in this case. Also, comparisons with the asymptotics of the exact 
closed- form expression, given in Eq. ( pil|) , are possible. 

The vanishing of the magnetic quantum numbers implies that (3, defined in Eq. ( p.l2[ ), 
is real. The case /3 = is simple because one of the j quantum numbers is then equal to 
the sum of the other two, and the integral in Eq. ( ^.111 ) may be evaluated exactly with a 
small amount of effort. Thus, we will consider the case (3 > 0. Two other facts that will be 
used throughout this appendix are that the set of quantum numbers is in the allowed region 
(since f3 is real) and that the j quantum numbers are integers (since ji — rrii is always an 
integer) . 

First, we work out the simplifications that occur in the all-real expression in Eq. ( p. 281 ). 
Equation (|3.24|) becomes 



TT / . . . 3 



(Bl) 



and Eq. ( p. 2^ ) becomes 



(Ji0j20|j0)^2 




cos 



TT 

L2 



Ui+J2-j) 



(B2) 
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This agrees with the first-order approximation of the exact expression, which can be obtained 
from the higher-order approximation [Eq. (P12|) ] of the exact result, given in Eq. (|B11|) . If 
Ji + J2 — j is odd, then both of the expressions are zero. (In this case ji + J2 + j is also odd 
since 2j is even.) If ji -|- J2 — j is even, then both have a sign of (— l)(Ji+-'2^-')/2. 

Next, we move on to the higher-order approximation. From Eq. ( |3.11| ) it is apparent 
that the two solutions for the cotangents of 6 and (p are related by simply reversing the 
signs. It follows from Eqs. ( |3.14D , ( [4.11| ) and ( [4. 121 ) that the values of the second and fourth- 



order derivatives of g{9, (p) are unchanged, while the third-order derivatives have their signs 
fiipped. Equations ( |4.7| ) and ( [4.1UD imply that the quantities ^4 and 5^ are unchanged. 



Finally, Eq. ( |3.17| ) implies that the quantity e^'^^°''^°'^ gets multiplied by (— l)Ji+-?2-i for the 
second root. Because j is an integer, this phase factor is the same as (— 1)-'i+j'2+j\ Since we 
are in the allowed region, we must sum over both stationary-phase points. We see that for 
odd values of ji + J2 + 3 the result is zero, while for even values of ji -|- J2 + 3 the result is 
twice the contribution obtained from one of the stationary-phase points. 

The values of the second derivatives of g{0, 0) at the stationary-phase points are obtained 
from Eqs. ( |3.11| ) and ( p.l4[ ), and the results simplify quite a bit. 



4j2 (ji+j) ^^^. 



ji + 32 + 3 



31+32+3 

From these equations results the following expression for the determinant of the 2x2 Hessian 
matrix of second derivatives of g{9, (j)). 

A ^ ^^3 I6J1J2J ,^„-, 

det ^,„ ,,„ = : : (B6) 

d{e,(f)y 31+32 + 3 

We note that this result is nonzero, and it does not vanish in any special cases that have 
nonzero j quantum numbers, which seems to contradict the statement made at the end of 
Sec. fTT] about the vanishing of the determinant when /3 vanishes. The resolution of this 
apparent contradiction has to do with the fact that in the stationary-phase analysis of this 
paper we do not simultaneously consider the cases mi = 1712 = m = and (3 = 0. These 
two conditions together would imply a is zero. The quantity a is defined in Eq. ( p.25|) . It 



vanishes when j = ji + J2, ji = 32 +3 or j2 = 3 + ji- In general, as long as a is nonzero, the 
expression for the determinant can be put (after some work) in the form 

d^g ^ (3P^+(3'P2 
^"'9(^,0)2 a^U +31+32?' ^ ' 

where Pi and P2 are (large) polynomials in the quantum numbers. This equation justifies 
the statement that the determinant is zero in cases where (3 is zero. On the other hand, in 
cases where mi = m2 = m = 0, Pi vanishes and a and f3 are equal, and the result simplifies 
to that shown in Eq. (p6|). The case of mi = m,2 = m = and /? = requires a separate 
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treatment. It is necessary to go back to the original integral representation for the Clebsch- 
Gordan coefficient. The integral may be approximated by the methods of stationary phase, 
but it is simpler just to evaluate it or Eq. ( |2.15| ) exactly, which is possible at that point. 

Higher-order derivatives of g{9, 0) at the stationary-phase points simplify as well. Equa- 
tions ( |4.7| ) and ( ^.10| ) for 84, and 5^ yield results that are much simpler than for the general 
case of nonzero magnetic quantum numbers. 



h + k = (jl^ J2 - 2 Ji=^ J2^ + ]1 ]2^ + ]1^] - Jl^ J2' J - Jl' J2^ J + J2^ J - Jl^ J2 f 



Wjl'j2'f 



3 -3 



2 aS 



3 -3 



- 3132' 3' - 23^3 -31 32 3 -3132^3' - 232' 3'+ 313" +323l/il'23 313213') (B8) 



This expression may be rewritten in a more compact form, as explained after Eq. (|B12D . 
As discussed above, for odd values of j + ji + J2 the higher-order approximation of the 
Clebsch-Gordan coefficient vanishes identically. For even values of j + ji + J2, the result 
reduces to 



. Ji+J2-i Zl -\- i 



3 + 3i + 32 
2vr/3 ]J 3+31+32 + 1 
1 1 



;i + <^4 + Sq 

1 



1 + 



1 

24 Vj 
1 



2 2 

ji 32 3 + 31 + 32 -3 + 31 + 32 3 - 3i + 32 3 + 3i - 32 



, (B9) 



where we have approximated the factorials in Nj^rn^_j2m2jm, using the form of Stirling's 
approximation that is appropriate for this order. 



x\ 



72^x^6-^ I 1 + 



12a; 



(BIO) 



The exact value of the Clebsch-Gordan coefficient is (Ref. [0, p. 87; Ref. [plj ) 

3 + 3i + 32 odd 

/ (-J1+J2+J)KJ1-J2+J)'(J1+J2-J)' 

V (Ji+i2+i+i)' ■ I ■ I ■ 

/ -Ji+j2+j 'ji ( h-32+i \\ ( h+32-3 \\i( h+i2+i \\ 3 + 3i + 32 cveu 

(Bll) 



(Ji0j20|j0) 






and this may be approximated using Eq. ( P10|) . The result is 
(Ji0j20|j0) 




. 31+32-3 
-1 2 



22+1 
2nf3 



j+ji+J2 
j+jl+J2 + l 



1 _ JJ1J2 

^ /32 



j + ji + 32 odd 
3 + 3i + 32 even 



(B12) 



To the order that we are working, the higher-order stationary-phase result, given in 
Eq. ( P9|) , and the corresponding approximation of the exact result, given in Eq. (|B12|) , agree. 
Equation ( |B9| ) contains two factors that have the form of unity plus a small correction. If 
these are multiplied out and only the ffist-order terms are kept, the result is the factor of 

(l — ^^M^\ in Eq. ( pi2|) . This shows that Eqs. ( p9D and (|B12| ) are equivalent, and it also 

provides an alternative way of writing the expression in Eq. 
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forbidden subregion 

I 

II 

III 

IV 

V 
VI 



Table I 






choice of root 


sign function 


lai 


lower 


1 


A3 


upper 


C—iyi-mi 


A2 


lower 


(_l)jl-i+m2 


Ai 


upper 


} — l\Jl+J2-j 


A3 


lower 


} — l\j2-j-mi 


A2 


upper 


( — \'\J2+m2 


Ai 



largest A 



Table I: For each forbidden subregion, the choice of root in Eq. ( ^.11[ ), the sign function 
as in Eq. (3.30) and the largest A [which determines the form of x, as in Eq. (3.31)] are 



given. 
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Figure Captions. 

Fig. 1: An example of a choice of three j- vectors, demonstrating that a set of quantum 
numbers is classically allowed. 

Fig. 2: The triangle-allowed region and the classically allowed region, shown in the rrii- 
1712 plane, for the case of j-values in the ratio j : j'l : j2 = 4 : 2 : 3. The six forbidden 
subregions are labeled with Roman numerals. 
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EigJ. 

An example of a choice of three j -vectors, demonstrating that a set 
of quantum numbers is classically allowed. 




The triangle-allowed region and the classically allowed region, shown in 
the nij-mj plane, for the case of j-values in the ratio j : j, : jj = 4:2:3. 
The six forbidden subregions are labeled with Roman numerals. 
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